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Abstract: We analyze power-suppressed contributions to singlet pseu- 
doscalar 77 and 77' meson transition form factors. These corrections stem 
from endpoint singularities and help improve the agreement between QCD 
theory and the experimental data, in particular, at low- momentum trans- 
fers. Using the CLEO data, we extract information on the profile of the r]i 
and 778 distribution amplitudes in the SU(2i)f octet-singlet basis employ- 
ing both the one-angle and the two-angle mixing schemes. In the former 
scheme, we find good agreement with the CLEO data, while in the second 
case, our approach requires non-asymptotic profiles for these mesons. 



1 Introduction 

Control over power-behaved corrections in QCD processes is crucial for the 
correct interpretation of high-precision experiments in which intact hadrons 
appear in the initial and/or final states. Prominent examples are meson- 
photon transition form factors, as measm'ed by the CLEO collaboration ^ 
for the pion and the r] and and the recent JLab high-precision data for 
the pion's electromagnetic form factor j^- 

Because, theoretically, the dynamics of such exclusive processes involves 
the corresponding meson distribution amplitudes, one can extract crucial 
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information about the nonperturbative partonic structure of pseudoscalar 
mesons. In contrast to the hard-scattering amplitude, that can be system- 
atically computed within perturbative QCD and is specific for each process, 
hadron distribution amplitudes are universal quantities that encode the par- 
tonic structure of hadrons. Their computation requires the application of 
nonperturbative methods, like QCD sum rules with nonlocal condensates — 
introduced in [SlI^El and recently improved in [H] — to derive a realistic pion 
distribution amplitude complying with the CLEO data on the pion-photon 
transition at the la level [Jj. In addition, this type of pion distribution 
amplitude was recently jH] used in conjunction with fixed-order ^ and 
resummed |1HI12) Analytic Perturbation Theory to calculate the pion's elec- 
tromagnetic form factor providing very good agreement with the existing 
data. 

Alternatively, one can use the factorization QCD approach in order to 
extract the shape of the pion (pseudoscalar meson) distribution amplitude 
directly from the data. The reliability of the latter possibility depends, how- 
ever, on the way one deals with fixed-order perturbative calculations. To 
minimize the influence of (disregarded) higher-order contributions, while ap- 
proaching the kinematic endpoint regions of the process in question (where 
the nonperturbative dynamics dominate), one may use Borel resummation 
techniques and, by this way, estimate power-behaved corrections. Indeed, 
this type of approach was recently used to compute the pion-photon 
transition form factor and determine the pion distribution amplitude with 
results for the latter quantity close to the profiles determined with the non- 
local QCD sum rules just mentioned. 

In the present exposition, we apply this type of approach (called the RC 
method) to extract the distribution amplitudes of the 77 and r]' mesons with 
particular focus being placed on their gluonic content jj^. More details of 
our approach can be found in ^3 El E[ El E] • 

2 77-77' mixing schemes 

Physical pseudoscalar mesons tj, -q' are admixtures of SU{3)p octet (r/g) and 
singlet (r/i) states: 

/ rj \ ^ ( cosOp -sinOp \ ( Vs \ (y. 
\ri' J V sin^p cos^p J\m J ' ^ ' 

where dp is the pseudoscalar mixing angle in the octet-singlet scheme (for a 
review and further references, see |2()j). 
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On the parton level, the states 773 and r/i are given by 

r?8 A / sin6'/ -cosOi \ ( -^{uu + dd)\ 
m J V cos^/ sin^, J y ss J ' ^ ' 

with 6 J being the ideal mixing angle. 

Then, in turn, rj and r/' are admixtures of qq pairs (in a quark-flavor 
basis) expressed via 



V 

rj I \ smap cosa^ I \ ss 



[uu + dd) 



(3) 



where Op = 6p — 9i + tt /2 denotes the deviation of the mixing angle from 
the ideal one due to the Ua{^) anomaly — in contrast to the vector meson 
(j) — oj system with a„ ~ 0. Note that the flavor-singlet pseudoscalar state 
contains also a gluon component: "gluonium". To accommodate the gluonic 
component, one ha s to extend the mixing scheme to a 3 x 3 matrix with three 
mixing angles; i.e..j21j 

(rj\ /cos 6p cos 7 + sin 6^^ cos sin 7 — sin cos 7 + cos cos sin 7 sin(;/>sin7\ 
r;' I = cos 9 J, sin 7 + sin 6p cos (j) sin ^p sin 7 + cos 6p cos cos 7 sin cos 7 1 
t / \ — sin 6p sin (p — cos 6p sin (p cos J 

(4) 

where i is a Glueball state and G = \gg) denotes gluonium. Note that \rf) ~ 
\r]^) (because ~ mg due to the Gell-Mann-Okubo mass formula), so that 
the 1 771) admixture is small with practically no room for a \G) contribution. 
Hence, 7 = 0, and, as a result. 







( cos Qy 




\-\ 


sin 6'p cos q 






\— sin 0p sin 



smf 




(5) 



A physical state is then a superposition of the sort 

It/;) =x|Q)+y|5) + z|G), + + ^2 ^ ^ ^g) 

with components (in the quark-flavor basis) given by 

\Q) = -^{uu + dd) 
\S) = \s-s) 

\G) = \g9). (7) 
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The physical states rj and 77' are 

\v)=Xr,\Q)+yrj\S), W)=Xrj'\Q)+yr,'\S)+Zr,,\G) (8) 

with mixing coefficients 

+ 2/^ = 1 > ^1'+ yl' + = 1 (9) 

related to the mixing angles 

= cosap , = — sina^ (10) 

and 

Xr,' = COS (f> sin ttp , Uri' = COS (f) COS Up , Zrj' = sin (p . (11) 
The SU{3)f octet-singlet basis is provided by 

|?7i) = —^\uu + dd + ss) 
V3 

\ri8) = \uu + dd - 2ss) . (12) 

v6 

One notes that the octct-singlet and the quark-flavor basis arc equivalent, 
but that the parameterizations of the decay constants are different. 

Let us now have a closer look to the decay constants of M = 77, 77' mesons. 
Their parameterization is defined via 

{Q\Jl,\M)=irpp^, (13) 

where J^g is the axial- vector current {i = Q, S or i = 1,8). In the quark- 
flavor basis, the decay constants follow the pattern of state mixing, i.e., 

fv = fQ fv =-fs sin Op 

= /q sin ap = fs cos ap . (14) 

In the octet-singlet basis the situation is different: 

= fsCOsOg /^ = -/isin6'i 

4, = /gsin^s ^^ = /iCOS^i. (15) 
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Note that in general, Q\ Op. In the present analysis, we use the 

octet-singlet basis with the one-angle (standard) parameterization with 

U = O.lSlGeV /i = 1.17/^ 

fs = 1.26/, Op = -15.4°, (16) 
and the two-angle mixing scheme with the parameters 

Op = -15.4°, Oi = -9.2°, ^8 = -21.2°. (17) 

3 Electromagnetic 777, 77^7 transition form factor 

In the Standard Hard-scattering Approach (HSA), the transverse momenta 
are neglected (collinear approximation) and the meson (vr, rj, -q' . . .) consists 
in leading twist (t = 2) only of valence \qq) and \gg) Fock states. Let us 
summarize some important issues: 

• The rj — rj' system shows flavor mixing due to the SU{3)p symmetry 
breaking and the U{1)\ axial anomaly. 

• The quark-singlet \r]i) and the gluonium state \gg) mix under evolu- 
tion; both carry flavor-singlet quantum numbers. 

• The gluon content of the rj' can reach the level of 26% |21j . 

• The meson-photon transition form factor contains a singlet and an 
octet part: Fu-yiQ^) = F^^{Q^) + F^^{Q^). 

• The singlet part, F^^{Q'^), has a quark and a gluonic component 
(HH — using the standard HSA; jTl] — using the endpoint-sensitive RC 
method). This means 

FM',iQ') = iTkM'^f^l) Tl{x,Q^^^l)) ( ^'l^'^^l ] , (18) 

where <^ denotes integration over longitudinal momentum fractions x 
from to 1 and /xp is the factorization scale. 

The transition form factor in the HSA can be expressed in terms of the 
convolution 

Fu'yiQ'^) = dx(/)M(2;,7i|)rH(x,Q^;/x|,/i|) (19) 
Jo 
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with = —q^ > and q is the four-momentum of the virtual photon. 
Figure ^ shows an example of the Feynman diagrams contributing to -FM7 
at NLO. 




Figure 1: A sample of a Feynman dia- 
gram contributing to the (pseudoscalar) 
meson-photon transition form factor at 
NLO. 



In Fig. ^ the partonic subprocess 7 * -I-7 q + q \s described by the 
hard-scattering amplitude rH(2;, Q^; /x^) (/iR being the renormalization 
scale), whereas the nonperturbative dynamics is contained in the universal 
meson distribution amplitude (/)M(a;,//|). Note that c^s satisfies a scalar 
evolution equation, analog to the vr case, while (j)i and (pg evolve together 
via a (2x2)-matrix evolution equation. Thus, we have |23j 



d(j)8{x,IJ,l) 

din fj.p 



(20) 



with a LO solution given in terms of Gegenbauer polynomials (x = 1 — x): 



h{x,^il) = 6x(l 



n=2,4,. 



(211 



and 



mx,fj-F) 



In Eq. (|2H) . the projection coefficients encode the nonperturbative in- 
formation that is not amenable to QCD perturbation theory, as we have 
already mentioned. On the other hand, the singlet (pi and gluonium (pg 
distribution amplitudes fulfill the matrix evolution equation 



(x,M,as(/x|)) ® 



^(pi{u,^l)^ 
^cpg{u,fil)^ 



(22) 
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with LO solutions provided by 



4>i{x,fj,p) = 6x(l - x) 



1+ j;i?M)C^/'(2x-l) 



n=2,4,. 



,(x,^|) = x(l -x)Y^ 5^(4) C^/_'i(2x - 1) 



n=2,4,... 

The normahzation conditions are 

/ dx(pi^s{x,fJ'F) = 1 
Jo JO 

The quark component of the singlet state reads 



{x,f4) = 0. 



(23) 



(24) 



01 (a::, fip) — 6xx 



1+ E 



(25) 

with the symmetry condition i;^i(x,^p) = 0i(a;,/i|), whereas the gluon 
component is 



Qs(mo) 



(26) 



with the symmetry condition cl)g{x,fip) 
anomalous dimensions jl7j are 



i{x,n^). The associated 



7S = 5 



T,", + r,, ± \/(7?,-7;,P + 47S,7?,J , (27) 



with 



7, 



99 



.7+ -7, 



99 



7, 



rn 



1 7: 



99 



91 



3 + 



6 7™ — 7, 



99 



n+1 



(n + l)(n + 2) —J 



(28) 



7, 



99 



n+l 



(n + l)(n + 2) 



(29) 



7 



^ 12n/ ^ + 3) . . 

(n + l)(n + 2)' ^s'i ■^3(n + l)(n + 2) ■ ^ ' 

The numerical values of these parameters (n/ = 3) are 



2 _ 2 2 _ ^ 2 _q 

Iqq — g ' Tgc; — -LJ-) 7c,g ~ 27 ' ~ 

2 48 2 101 , 16 3 1 

7+ - -y, 7---^, P2-y, ^2--^- (31) 



The required Gegenbauer polynomials are 

Cl^'^{x-x) = - [5(x-x)2-l] =6(l-5xx) 

Cl''^{x-x) = ^{x-x). (32) 

4 Hard-scattering amplitudes for the 777 and 77^7 
transition 

The form factor for the 777 and r]'^ transition, given by Fm7(Q^) = -^M'y(Q^) + 
F]y[^((5^), contains a singlet part comprising quark and gluon components: 

[^H,NLo(*'<3^MF)<^^'l(a;,/UF) +^^,NLo(2^'<3^MF)0s(a;,MF)]} ■ (33) 

The octet part contains only a quark component; it reads 



Jo 



^H,Lo(^)<^8(a;,MF) 



The expressions for the involved hard-scattering amplitudes are 

-1 , ™-i 



(34) 



r^,Lo(^)=^ ; (35) 



-^H,NLo(^' Mf) 



2 X In X 
In a; 



1 02 
+ -(3 + 21nx)ln-^ + (a; 

x fj,^ 



(36) 
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^H,NLo(^' Q ' 4) = + (^6 - -J Inx + In ^ - (a; ^ x) (37) 

and the charge factors read 

Ni = ^{el + el + e',), Ns = ^ {el + ej - 2e^,) . (38) 

5 ?77, 77^7 transition form factor in the RC approach 

Let us outhne here the essentials of the endpoint-sensitive RC method. 

• Solve the renormalization group equation for as{XQ^) in terms of 
"s(Q^) UM to Q^s(Q^) accuracy. 

• Expand the hard-scattering amplitude T{Q'^) of the process as a power 
series in as{Q'^) with factorially growing coefficients Cn ^ {n — 1)!. 

• Use the Borel integral technique to resum them by 

— determining first the Borel transform B[T] of this series 

— inverting then B[T] to get 

/•oo 

Jj.]resum^g2) ~ RV. / duexp 

Jo 



-Attu 



/3oas(Q2) 



B[T]iu) 



At this point a couple of important remarks are in order, (i) The 
Borel transforms contain poles on the positive u axis that are exactly 
IR renormalon poles; hence a principal value (R-V.) prescription has 
to be used, (ii) A direct way to obtain the Borel resummed expressions 
is via the Inverse Laplace Transformation. 

Then, one finds 

poo 

a^(xQ^)= due-''^R{u,t)x-'' (39) 

Po Jo 

with 

i?(u,t) = l-^n(l-7E-lnt-lnn). (40) 
Po 

• Endpoint singularities x — > 1 x ^ 1 transform into IR renormalon 
(multi-)pole divergences at uq = n (in our case n = 1,2,3,4) in the 
Borel u plane. 
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• Removing these poles via the P.V. prescription, we obtain resummed 
expressions for 



• The pole at uq = n of the Borel plane corresponds to power-suppressed 
corrections ~ (l/Q^)" contained in the scaled form factors. 

Let us close this section, by commenting upon the importance of power 
corrections from the theoretical point of view and in comparison with the 
standard HSA. The latter prefers to set fi'^ = Q^. Then, large NLO loga- 
rithms are present. The RC method sets instead ju!^ = xQ^. As a result, the 
term \n{iJ,^/xQ'^) in the NLO contribution is eliminated, but the integration 
over X gives rise to power-suppressed contributions in the endpoint regions 
X ^ 0,1. Note in this context that because asymptotically both approaches 
have to yield the same results, one has to verify that the induced power 
corrections do not affect this regime, leaving the asymptotic behavior of 
perturbative QCD unchanged. Hence, in technical terms, one has to ensure 
that 

/ due-''^R{u,t)^^^ due-'^K 
Jo Jo 

In view of the above remarks, the best (perturbative) procedure is the one 
that minimizes the NLO contribution while keeping power corrections under 
control. 

The present analysis employs the following scales: 

• /x^ = xQ^ (renormalization scale) 

• a(»/=4) = 0.25 GeV 

• //q = 1 GeV^ (normalization scale) 

• iJ,p = Q'^ (factorization scale) 

The estimated influence of higher-twist uncertainties is of the order of (10 — 
15)%. 
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6 Borel Resummed 777 and 77^7 transition form fac- 
tors 



The NLO expression for the transition form factor, calculated with the RC 
method comprises a quark component 



g2F]li^(Q2)quark ^ a,{Q'' x)t{x , fji'p) , I^'f) 

+ as{Q'^x)t{x,iJ,p) (f>i{x,fj,p) 
= 2as(Q^x)t(x,/i|) 0i(a;,|ip) 



with 



t{x,np) 



1 



,9 x In X 
In^x 9 



1 

+ -(3 + 21nx)In-^ 
X nj, 



and a gluon component 

M7 



with 



l{x,ii^) 



2^ xln'^x 4:\ ^ ^ fQ^\ xlnx 



+ 6 



In x + 2 In 



(41) 
(42) 

(43) 
(44) 



Summing up, we can write-in the context of the RC method-the transition 
form factors Q'^F^^{Q'^) and Q'^F^^{Q'^) as follows 



< 01 (x, ^f) + as{Q'^x)g{x, fip) (j)g{x, f/p)] 



(45) 



and 



fuNs 



Tj^ g(x) (g) (/>8(x, /x|) + ^ as((3^x)t(x, /i|) 



(46) 



Now recall that the running coupling as{xQ ) in terms of as{Q ) [23 reads 



asiQ'^x) 



1 + Inx/t 



as{Q^)f3i ln[l + lnx/t] 
27r/3o 1 + lnx/t 
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where we have used 
in 



q2 2 19 

2\ -^^~\2^ /3o = 11 - i^nj , /3i = 51 - —Uf . 



In this way, we finahy arrive at the following expressions for the re- 
summed singlet and octet transition form factors within the RC method: 



:i+^(mI)) 



/"CXD /"CXD 

X / due-''*R{u,t)Qi{u) - 5A{fil) / due'''* R{u,t)Q2[u 
Jo Jo 

+ /"dne-"*i?(u,t)G(tx)l , ( 

Po Jo J 



PO 



/■OO 

'l + C{^ll)) / due- 
Jo 



■ut 



f'OO 

X R{u,t)Qi{u)- 5C{nl) / due-''^R{u,t)Q2{u) 

Jo 



■ (48) 



In the above expressions, the following abbreviations have been used |14j : 

R{u,t) = l-'^u{l-jE-l^t-lnu), (49) 
Po 



Qi(n) 



2a 1 + 2a a - 3 



(l-u)3 (2-tx)3 (l-^/)2 ' (2-u)2 ' "(l-n)(2-^/) ' 



Q2(n) 



4 2 

+ 



2a 



+ 



(2-n)3 (3-n)3 (4-n)3 (2 - n)2 

l + 4a l + 2a a-3 
+ 6 



(3-n)2 (4-n)2 (2 - u)(3 - u)(4 - n) 



(50) 



G(n) 



2 2 

+ 



5 — a , 3 — a 
2— 



(4-n)3 (3-n)3 (2 - u)2 (3 - u)2 (4 - n) 



(51) 
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Figure 2: 77 (lower curves) and 77' (upper 
curves) transition form factors obtained 
with the standard HSA in The 
shaded area corresponds to the range of 
values Bf(^§) = -0.04 ±0.04 B^{i4) = 
-0.08±0.04 Bf (a^o) = 9±12 (see [22123 
for further details). 



with a = ln{Q / Hp) and the Gegenbauer coefficients being given by 



6BI 





48 

81 






101 

81 


_as(Ato). 




15 


.as(/io). 





Bill 



16BI 



48 
81 



+ 5BI 



101 

81 



6Bl 



Qs(^l) 



50 
81 



(52) 



7 Phenomenological Analysis 



In this section we perform numerical computations of the Borel resummed 
and rescaled by 777 and "q'^y transition form factors in order to extract the 
T] and 77' meson distribution amplitudes from the CLEO data. We shall also 
compare our theoretical predictions with those obtained with the standard 
HSA (221 ESI I being to reveal the role of power corrections at low- 

momentum transfer in the exclusive process under consideration. 

Let us start our discussion by quoting the results obtained in |22] (see 
also |26j ) using the standard HSA. Their main predictions are shown in Fig. 
[21 in comparison with the CLEO data pQ. 

As one sees from this figure, the agreement between the theoretical pre- 
dictions and the low-momentum data is rather poor — especially when us- 
ing asymptotic profiles for the 77, 77' meson distribution amplitudes. To 
decrease the magnitude of the form factors at low Q^, and achieve this 
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r\y Bj'(ti,)=b/(ii,)=o.o2 




w 0.3- 

"a 


Ti'y B;(,i,)=B;(n,)=0.02 








0,2- 










0.1- 




(b) 



5 10 15 20 5 10 15 20 

Q' (GeV') Q' (GeV") 



Figure 3: Predictions for the scaled form factors as functions of of the 777 (left 
panel) and rj'^ (right panel) electromagnetic transition. For the solid curves the 
designation is -B|(?7i) = 0. The dashed lines correspond to i3f (r^i) = 10; for the 
dash-dotted curves we use -Bf ('yi) = 15. The data are taken from Rcf. (J. 



way a better agreement with the data, the standard HSA would call for 
the two-angles mixing scheme and for distribution amplitudes mainly with 
-B2(??i), B2{ris) < 0. The inclusion of power-law corrections changes the low- 
momentum behavior of the form-factor predictions significantly, as one ob- 
serves from Fig. |3| Indeed, using the standard octet-singlet mixing scheme, 
one can reproduce the trend of the CLEO data rather well in the whole 
momentum range explored — especially with a non negligible gluon contri- 
bution (the Gegenbauer coefficients are given in Fig. ^ — because the effect 
of power corrections is to enhance the absolute value of the NLO correction 
to the form factors by more than a factor of 2.5 — 3. Since the contribution 
of the NLO term to the form factors is negative, the power corrections re- 
duce the leading-order prediction for the form factors considerably, while at 
the highest values measured by the CLEO collaboration this influence 
becomes more moderate. 



The la regions in the form of shaded areas for the scaled form factors 
for the rj'f and ri'^y transition in the RC method and using the octet-singlet 
scheme are displayed in Fig. \^ The central line corresponds to the coeffi- 
cients values 5|(r?i) = B^{r]s) = 0.05; 5f(??i) = 17. A full-fledged dis- 
cussion of these issues is given in ^1], together with error estimates arising 
from varying the values of the theoretical parameters used in the analysis. 
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(a) 



(b) 



in 



2(1 



15 



20 



q' (GeV') 



(GeV') 



Figure 4: The 777 (left) and (right) scaled transition form factors as functions 
of Q^. The central solid curves are found using the values B\{r\]) = B\(t]^) = 0.05 
and B?^[r]\) = 17. The shaded areas demonstrate the Icr regions for the transition 
form factors. 



It is important to emphasize that our calculations do not exclude the 
usage of the two-angles mixing scheme in conjunction with the RC method. 
But in such considerably larger contribution of the non-asymptotic 

terms to the distribution amplitudes of the r\\ and r^g states would be re- 
quired. Carrying out such a computation |14j . we obtained the results shown 
in Fig. 121 Inspection of the left panel of this figure reveals that the r/7 tran- 
sition FF found within this scheme lies significantly lower than the data. 
Therefore, to improve the agreement with the experimental data, a rela- 
tively large contribution of the first Gegenbauer polynomial to the distribu- 
tion amplitudes of the r]\ and r/s states seems necessary. The Gegenbauer 
coefficients corresponding to the predictions shown in Fig. [S] are B\ (rji ) = 
0.15, S|(r/8) = 0.15 and i?2(??i) = 18. We consider the values i?2(^i) = 
i?2(??8) = 0.15 as actually determining the lower bound for the admissi- 
ble set of distribution amplitudes in the context of the two-angles mixing 
parameterization scheme. Hence, in that scheme, we obtain 



Blim) = Blirjs) = 0.15, Blim) G [16, 20] . 



(53) 



15 



s 



ny 






















(a) 





5 


10 


15 2( 













n 

(b) 



Q (GeV ) 



Figure 5: The 777 (a) and rl^ (b) electromagnetic transition form factors vs. . 
The solid lines correspond to the ordinary octet-singlet mixing scheme with param- 
eters B\(r\\) = SKryg) = 0.02 and i3|(77i) — 18. The broken lines are obtained 
within the two-angles mixing scheme. The dashed lines describe the situation with 
the same parameters as the solid curves. The parameters for the dash-dotted curves 
are (r,i) = BUrjs) = 0.15, B|(?7i) = 18. 



Let us close this section by summarizing the main differences between the 
standard HSA and the RC method: (i) Form factors in the HSA overshoot 
the CLEO data — especially in the low region — even with the NLO correc- 
tions included, (ii) Values of the Gegenbauer coefficients -B|(r/i), i?2(^8) > 
increase the disagreement, while 51(^1) > reduces the disagreement. 
Hence, a better agreement with the CLEO data would call for the two- 
angles mixing scheme and B2(??i) ,S2(%) < 0- (iii) The inclusion of power 
corrections enhances the (negative) NLO correction to the form factors at 
low by factors 2.5 — 3. In order to quantify these statements, we show 
in Fig. El the numerical results for the ratio 



IHSA 



[Q2Fm^(Q2)] 



(54) 



NLO 



for some selected values of the expansion coefficients. As a result, the RC 
method, employing the one-angle mixing scheme, is in good agreement with 
the CLEO data, (iv) Using instead the two-angles mixing scheme, the RC 
method favors non-asymptotic profiles for the distribution amplitudes of rji 
and r/8, e.g., i?2('?i)) ^2(^8) ^ 0-15 and Sfl^i) ^ [16, 20], while the region 
-6|(Q^ =^ 2 GeV^) < seems to be incompatible with the CLEO data. 
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Figure 6: The ratio -R(Q^) for the 

?77 form factor. The solid line corre- 
sponds to the input parameters B2{r]i) = 
B^irjs) = Bf(r/i) = 0. The dash-dotted 
curve describes the same ratio, but for 
Blivi) = B^im) = 0, Blini) = 14, while 
the dashed one corresponds to ^2(771) = 
B|(7?8) = 0.05, Bf (771) = 10. 



8 Conclusions 

The renormalon-inspired RC method enables the inclusion of power correc- 
tions originating from the kinematic cndpoint region (.x 0, 1), where non- 
perturbative QCD dominates and fixed-order perturbative computations of 
such corrections yields divergent results. We found that power-suppressed 
ambiguities to form factors vary between 3% at high and 11% at low 
values. On the other hand, wc have verified that the asymptotic limit of 
[Q'^FM'y{Q'^)Y'^^^'^ coincides, as it should, with the standard HSA result, 
leaving the asymptotic properties of QCD perturbation theory unchanged. 
The effect of power corrections at < 5 GeV^ enhances the (negative) NLO 
correction by 2.5 — 3 times, providing this way agreement with the trend of 
the CLEO data. In the standard octct-singlct scheme we found i3|(r/i) = 
-^ll^) — 0-055 lb 0.065, -B|(r/i) = 18 it 4.5, whereas in the two-angles 
mixing scheme, we found = B^{rfs) > 0.15, 5f (r?i) e [16, 20]. The 

distribution amplitude of the r] and rj' mesons, obtained in this work, can be 
useful in the investigation of other exclusive processes that involve rj and rj' 
mesons, especially at lower momentum-transfer values, where the standard 
HSA is most unreliable. 
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